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Abstract 

We discuss T-odd correlations between jet and lepton momenta in ti events at the Tevatron that 
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CsJ can be used to search for CP violation. We identify correlations suitable for the lepton plus jets 
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the top-quark anomalous couplings, including its chromo-electric dipole moment, and we estimate 
the limits that can be placed at the Tevatron. 
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I. INTRODUCTION 



Looking for new sources of CP violation remains one of the important goals of high energy 
colliders. Processes that have been considered before include production and decay of top 
quark pairs [l]-[6j], production and decay of electroweak gauge bosons [zj, and production and 
decay of new particles jgj]. 

The Tevatron has now observed hundreds of ti events and will eventually collect a few 
thousand. With this in mind, it is interesting to revisit the issue of possible CP violation 
in these events. Under optimal conditions, a sample of a few thousand events would have a 
statistical sensitivity to a CP violating asymmetry at the few percent level. 

In this paper we study T-odd triple product correlations of the sort first discussed in 
Ref. 0, [9| and find observables suitable for ti events in which the top quark pairs decay into 
a lepton plus jets or purely hadronically. Of particular interest are observables that only 
require reconstruction of two b jets (but there is no need to distinguish between b and b); 
one or no hard leptons; and non-6 jets ordered by pt- 

It is well known that CP violation in the standard model (SM) is too small to induce a 
signal at an observable level in the processes we consider. We will discuss two CP violating 
effective interactions that should serve as benchmarks for the sensitivity of the Tevatron to 
CP violation in ti events. 

As in our LHC study |(3], we parametrize CP violation using anomalous top-quark cou- 
plings. The ti production process is modified relative to the SM by the chromo-electric 
dipole moment (CEDM) of the top-quark via the interaction 

C-cdm = -igs^ta^G^t, (1) 

where g s is the strong coupling constant and G^ u is the usual gluon field strength tensor. 
The CEDM is induced in principle by any theory that violates CP, and estimates for its 
size in several models can be found in Ref. Typical estimates presented in Ref. [l[ for 
the size of d suggest that it may be too small to yield observable signals at the Tevatron. 
Nevertheless, we view the study of this coupling at the Tevatron as a valuable preliminary 
to future LHC studies. 

We also consider CP violation in the decay vertices t — > bW + and i — > bW~~ via the 
anomalous coupling / defined by 1 , 

T» wth = --tLv*u(p b ) [^PL-ife^f+^a^ipt-p^P^uipt) 

f ™ = "7! Vtb §{pi) i lflPL ~ l f et{ ~ H+Sf)a " u (pt ~ n)uP L ] v( P - b ), (2) 

In Eq. [2] we have explicitly split the phase of / into a CP violating phase 0/ and a CP 
conserving, unitarity, phase 6f. 



1 As discussed in the literature, other anomalous couplings will not interfere with the SM and we will not 
consider them here 10 1 . 
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The CP violating anomalous couplings, d and / sin </>/,_ have been recently revisited vis-a- 
vis the upcoming LHC experiments in Ref. [5| and Ref. p]. In Ref. jij], general results were 
derived for the T-odd correlations induced by these two couplings for both gluon fusion and 
light qq annihilation ti production processes. In the appendix we specialize those general 
results to the specific processes that are relevant for the Tevatron. In Ref. jf! a numerical 
analysis was carried out for LHC concentrating on the dilepton signal, which is not viable 
at the Tevatron due to the small number of events. The new signals we discuss in this paper 
pertain to the lepton plus jets and all hadronic decay modes of the top-quark pairs and can 
also be used at LHC. 



II. OBSERVABLES 

In Ref. jij all the linearly independent T-odd correlations induced by anomalous top- 
quark couplings were identified. From these we need to project out the ones that are most 
suitable for the Tevatron and two considerations come into play. The first one, already 
discussed in our application to the LHC in Ref. [sj], is that we want to use only momenta 
that can be reconstructed experimentally. The second one is that, due to the low statistics 
at the Tevatron, we will be dealing with at least one hadronic decay of the W boson. 

We will consider the following correlations 2 : 

• For the lepton (muon) plus jets process pp — > ti — > bbfijij 2 + $t'- 

/n t \ tt CM ^ / ^ ^ \ 

Oi = e\Pt,VuPb,Pb) > oc p t ■ [p b x p- b ) 

e> 2 = e(P,p b +p h p e ,p jl ) — > oc (Pb + Pi) ■ (Pe x pji) 

03 = Qet(Pb,PhPe,Pji) bb CM > oc QePb- (pe x Pji) 

4 = Qe e(P,p b - pi,Pe,Pji) ^ oc Q e {pb - Pi) ■ (pe x Pji) 

Ol = q-{Pb~ Pb) Q, Pb, Pb) ^ OC Pbeam ' {Pb ~ Pb) Pbeam ' {Pb X p~ b ) ■ (3) 

• For the multi-jet process pp — >• ti — >■ &&J1J2J1' J2' : 

O x = e(p t ,pt,p b ,pi) U ) oc p t ■ {p b x p- b ) 

05 = £{Pb,PhPji,Pji') ► oc p b ■ (Pji x Pj V ) 

6 = e(p b ,p- b ,p jl +p j2 ,Pji> +Pjv) ~ ° M) oc [p jX + p j2 ) ■ (p b x pi) 

Ol = Q- {Pb- Pb) € i P , Q, Pb, Pb) ^ OC Pbeam " (Pb ~ Pb) Pbeam ' (Pb X Pb) ■ (4) 

In Eqs. [3] and H] we have shown two expressions for each of the correlations. The first 
one is valid in any frame and in particular can be used in the lab frame. The second one 

2 Here we use the Levi-Civita tensor contracted with four vectors e(a, b, c, d) = e llva pa fJ 'b 1 ' 'c a dfi with the sign 
convention £0123 = 1- We also use s,t,u to refer to the parton level Mandelstam variables for qq — > ti. 
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shows the correlation in a particular frame in which it reduces to a simple triple vector 
product. In these expressions P is the sum of the proton and antiproton four- momenta; 
q is the difference of the proton and antiproton four-momenta; p b i refers to the b or b jet 
momenta; pi refers to the momenta of a lepton that has been identified as originating from 
t or f decay in lepton plus jets events (in our analysis we only consider muons); pji, pji 
refer to non-6 jets ordered by p? (hardness) that reconstruct a W; primes denote the two 
jets associated with the second W in the all hadronic case. In 56 it is not necessary to 
distinguish b and b jets. It is only necessary to associate ji^ with one of the b jets and j[ 2 
with the other one when reconstructing the top-quark pair event. The hardness of the jet is 
defined in the usual way, the hardest jet being that with the largest transverse momentum, 
i.e. pTj > PTj 2 , and Qe is the lepton charge (for some of the monoleptonic signals, lepton 
charge id is needed). 

Notice that some of the correlations require differentiating between the b and b jets 
but others don't. In addition, 0\ requires the reconstruction of the top momenta. This 
correlation is the one closest to the form that appears in the parton level calculation, and 
in a perfect reconstruction situation it is identical to Oq. In addition, two of the examples 
given, 0\ and can be used for both processes. Finally we note that there are many other 
possibilities that we have not listed. 

All the correlations listed above are CP odd, as can be seen most easily in the specific 
reference frames given. For example, O2 in the lab frame becomes 

2 ^ VS [(p b + p- b ) ■ x p n ) + (p^- x p 3l ))] 

^ {-)VS [{p b + pi) ■ (Ov x p n ) + (p^ x p n ))] . (5) 

Eq. Oalso clarifies what is meant by O^'- events with a /i + and a W~ decaying to two jets 
will contribute to the first term in the sum in the first line. Events with a pT and a W + 

CP 

decaying to two jets contribute to the second term. The assignment pj\ — > —Pji on the 
second line of Eq. Estates that if CP is conserved, the probability for a given jet originating 
from a quark q in a two jet W + decay to be the hardest one, is equal to the probability 
for the corresponding jet originating from the anti-quark q in a two jet W~ decay to be 
the hardest one. These statements are verified in our numerical simulations both explicitly 
and by the fact that the asymmetry is induced by CP violating couplings but vanishes 
for CP conserving ones. In an experimental analysis it will be important to implement 
additional cuts in a way that is CP blind, typically requiring the same cuts for particles and 
anti-particles. 

Use of the lepton charge in some of the correlations allows us to construct CP odd and 
CP even correlations with the same set of momenta. We exploit this to construct the T-odd 
(but CP even) correlations sensitive to strong phases: 

O a = e(P,p b -pi,p £ ,p jl ) A oc (pb-Pi) ■ {pe x p n ) 
O b = Q e e(P, p b + pi, p e , Pji) oc Q e (p b + p- b ) ■ (pt x p^) 
O c = e(P,p b +Pb,Pji,Pji>) ^ oc {Pb + Ph) ■ (Pji x Pn')- (6) 
The first two have CP-odd analogues in C 4 and G 2 , respectively. 
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Our observables will be the lab frame distributions da/dOi for the correlations listed 
above, as well as their associated integrated counting asymmetries 

A = NeventsjOi > 0) ~ NeventsjOj < 0) (?) 

1 ~ N events (O t > 0) + N events (Oi < 0) ' U 

the denominator being just the total number of events in all cases. When our numerical 
results for the integrated asymmetries are very small we distinguish between very small 
asymmetries and vanishing asymmetries as described in Ref. jfij]. 



III. NUMERICAL ANALYSIS 

Our numerical study in this paper corresponds to the implementation of analytic results 
presented in Ref. j^] . The T-odd correlations for the parton level processes that are relevant 
at the Tevatron are not explicitly written in Ref. jB], so we present them in the Ap pen dix 
for convenience. The numerical studies are performed with the aid of Madgraph [ll|-|l3 
following the procedure outlined in Ref. For the lepton plus jets channel, we begin with 
the standard model processes qq (or gg) — > tt — > b£ + vbud (or bi~vbdu) implemented in Mad- 
graph according to the decay chain feature described in Ref. [13j|. This decay chain feature 
is chosen for consistency with the approximations in the analytical calculation of the CP 
violating interference term presented in Ref. jsj, in which the narrow width approximation 
is used for the intermediate top quark and W boson states. The expressions from Ref. ^ 



(Eqs. lAlHAllj) are then added to the spin and color averaged matrix element squared for the 



SM (which Madgraph calculates automatically) and the resulting code is used to generate 
events. A similar procedure is followed for the purely hadronic decay of tt with the relevant 
parton level processes. In this case both W's decay into a pair of quarks and we only con- 
sider the final states u, d, s, c without Cabibbo mixing. The code used to generate events 
is, therefore, missing the terms that are completely due to new physics: those proportional 
to the anomalous couplings squared. This approximation is justified because those terms 
do not generate T-odd correlations. In addition, as long as the conditions that allow us to 
write the new physics in terms of anomalous couplings remain valid, their contribution to 
the total cross-section is small. 

For event generation we require the top quark and W boson intermediate states to be 
within 15 widths of their mass shell, and two sets of cuts. The first set of cuts includes 
a minimum transverse momentum for all leptons and jets, a minimal separation between 
them, and a pseudorapidity acceptance range: 

p Tfl . > 20 GeV, p Tb i > 25 GeV, < 2.5, AR ik = - Vk ) 2 + (0, - fe ) 2 > 0.4 (8) 
with i,k — b, b,j, /i. 

For the second set of cuts (in the lepton plus jets channel) we add a missing transverse 
energy requirement 

$ r > 30 GeV. (9) 

We use SM parameter values as in Madgraph, except for m& = 0; and we use the CTEQ-6L1 
parton distribution functions. 
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A. Process pp — >■ tt — >■ hb[ij\22 + $t 

We first estimate the counting asymmetries by generating 10 6 events for each of the four 
cases: d = 5 x 1(T 3 GeV" 1 ; / sin0/ = 5 x KT 3 GeV" 1 ; / sin 5/ = 5 x 10" 3 GeV" 1 and 
d = f = 0. These cases correspond to CP violation in the production vertex, CP violation 
in the decay vertex, strong phases in the decay vertex and the lowest order SM respectively. 
The relatively large number 5 x 10~ 3 GeV -1 is chosen to facilitate distinguishing signals from 
statistical fluctuations. Once we establish a non-zero asymmetry we can cast our result as a 
function of the anomalous couplings since the asymmetries are linear in them. As mentioned 
above, we include the new physics only through its interference with the SM. Since these 
T-odd terms are also P odd, they do not affect the cross-sections as they integrate to zero. 
For this reason the total number of events is the same as in the standard model. 

The results with the set of cuts Eq. [8] are shown in Table fl] After these cuts are applied 
there remain approximately 6.7 x 10 5 events, leading to the 3a statistical sensitivity shown 
in the first column. The results show that all the CP-odd correlations vanish for the two 
CP conserving cases (SM and / sin 5/), and that the CP-even correlations vanish except for 
the CP conserving case with unitarity phases, / sin 5/. This establishes numerically that, 
at least at this level of sensitivity, there is no CP conserving contamination of the CP-odd 
signals or vice versa. 





3/y/N 


Ai 


A 2 


^3 


A 4 


A 7 


A a 


A b 


d 


3.7 


-66.9 


-37.4 


-100.6 


75.8 


40.4 


-3.4 


-1.8 


/sin <j)f 


3.7 


-7.2 


-60.8 


-8.2 


-36.7 


10.6 


0.6 


-1.9 


/sin (5/ 


3.7 


-0.8 


-1.0 





-1.7 


1.9 


-49.3 


-51.9 


SM 


3.7 


2.6 


-0.6 


0.4 


0.1 


0.2 


-0.4 


0.3 



TABLE I: Integrated asymmetries with cuts given in Eq. [8] for d, fsin(<pf,5f) 
= 5 x 1(T 3 GeV^ 1 in units of 10~ 3 , and the SM. 

In Table [TT] we show the same results with the additional missing Ej- requirement of 
Eq. [9j This additional cut further reduces the number of generated events to about 6 x 10 5 , 
and slightly decreases the statistical sensitivity. The effect of this cut is minimal on all 
asymmetries, making it very desirable for reducing background. 

Now we summarize our results for the asymmetries in the process pp — > tt — > bbfijij2 + $t 
with cuts given in Eqs. [HI E] in terms of the dimensionless anomalous couplings 

d t = dm t , ft = fm t (10) 



6 







M 


A 2 


^3 


A 4 


A 7 


K 


A 


d 


3.9 


-66.4 


-38.9 


-102.3 


76.5 


36.4 


-3.0 


-1.4 


/sin 4> f 


3.9 


-17.2 


-66.8 


-18.8 


-30.8 


7.0 


0.7 


-2.3 


/sin S f 


3.9 


0.4 


-1.1 


1.6 


-3.1 


1.0 


-44.1 


-56.8 


SM 


3.9 


2.5 


-0.5 


0.5 


-0.2 


0.4 


-0.3 


0.6 



TABLE II: Integrated asymmetries with cuts given in Eqs. El E] for d, /sin(0/,5/) 
= 5 x 1(T 3 GeV" 1 in units of 10~ 3 , and the SM. 

with m t = 171.2 GeV. We find, 

Ax = -0.078 d t - 0.020 / t sin 0/ 

A 2 = -0.045 dt- 0.078 / t sin 0/ 

A 3 = -0. 12 dt - 0.022 /i sin <p } 

A 4 = 0.089 dt- 0.036 / t sin 0/ 

A 7 = 0.043 d t + 0.008 /t sin 0/ 

A a = -0.052 f t sin 5/ 

A fe = -0.066 / t sin 5/. (11) 

In addition to the integrated counting asymmetries one can look for asymmetries in the 
distributions dajdOi. In Figured] we compare the distributions for da/d0 2 induced by (a) 
d = 5 x 10~ 3 GeV- 1 and (b) /sin0/ = 5 x 10~ 3 GeV -1 to the SM. 

It is instructive to discuss A 2 in some detail to understand the role of the hardest jet 
momenta. The lepton and (non-6) jet momenta that appear in this correlation act to some 
extent as the spin analyzers in the t and i decays. It is well known that the best spin 
analyzers in the top-quark rest frame are the charged lepton momentum (for semileptonic 
top decay) and the d-quark momentum (for hadronic decay) |14]. Of course, it is not 
possible to tag the d-quark jet in experiment, but at the event generator level we can see 
how things work. To this effect we define A 2 d, u , the counting asymmetry corresponding to 
2 d,u = t(P,Pb + Phi Vti Pd,u) ■ These asymmetries, as the original A 2 , are interpreted as the 
sum of processes with fi + , d( or u) from semileptonic decay of t and hadronic decay of i, and 
processes with d( or u) from semileptonic decay of t and hadronic decay of t. With the 
cuts in Eqs. [8] and [9] and with^d = 5 x 10~ 3 GeV~\ we find A 2d = 1.5% and A 2u = -2.7%. 
Interestingly, the asymmetry A 2 d, is sma//erthan A 2 = —3.9%, which appears in Table HT1 To 
understand what happens, we show in Figure [2] the differential distribution of the numerator 
of A 2 with respect to r, the ratio of d-quark transverse momentum to w-quark transverse 
momentum in t — > bud or t — > bud decay. As can be seen in the figure, if one chooses 
the d-quark momentum in the lab frame to construct this particular correlation, there is a 
partial cancellation between the regions with pxd > Ptu and pxu > PTd- This cancellation 
is removed by choosing instead the hardest jet resulting in the larger A 2 . The fact the A 2 
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SM 



•1 -0.5 
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FIG. 1: Differential distributions da/d0 2 for (a) d — 5 x 10 3 GeV 1 and (b) 
/sin0/ = 5 x 10~ 3 GeV -1 compared to the SM. 2 is given in units of m\. 



0.02 




FIG. 2: Differential distributions of s, the numerator of A 2 d (or A 2u ) , with respect to r, 
the ratio of d-quark transverse momentum to w-quark transverse momentum in t — >■ bud or 

t — »■ 6wci decay. 



is larger when using the hardest jet instead of the <i-quark jet appears to be unique to this 
correlation. 

Using our generated events, we estimate that the <i-quark jet has a larger pt than the 
u-quark jet 44% of the time. We can also verify that within statistical errors, the probability 



S 



of d being the hardest jet in t — > bud is indeed the same as the probability of d being in the 
hardest jet in i — > bud. 

B. Process pp — > tt — >■ bbji j^ji'h' 

In this case we only use the cuts of Eq. [S] as there is no missing energy. The results are 
shown in Table II III for about 7.4 x 10 5 generated events. As expected, the signals A% and 
Aq are the same at the parton level and it remains to be seen what dilution there is after 
hadronization. 





3/VN 


A l 


A 5 


A 6 


A 7 


A c 


d 


3.5 


-61.2 


-54.6 


-61.2 


38.8 


1.1 


/sin (p f 


3.5 


-7.1 


-7.8 


-7.1 


5.8 


-1.0 


/sin 5 f 


3.5 


-1.8 


-1.5 


-1.8 


-0.5 


9.6 


SM 


3.5 


0.7 


0.5 


0.7 


1.0 


1.1 



TABLE III: Integrated asymmetries for signal 2 with cuts given in Eq. [8] for d, /sin(0, 5) 
= 5 x 1(T 3 GeV^ 1 in units of 10" 3 , and the SM. 

Using our results In Table Hn] for the process pp — » ti — )■ bbj^jvjy we find 

A 1 = -0.072 ck - 0.008 f t sin 0/ 

A 5 = -0.064 dt - 0.009 f t sin (j) f 

A 6 = -0.072 dt - 0.008 / t sin <f> f 

A 7 = 0.045^ + 0.007^^0/ 

A c = 0.011 / t sin 5 f . (12) 

In Figure E] we compare the distributions for da/dO^ induced by (a) d = 5 x 10~ 3 GeV -1 
and (b) fsm<f> f = 5 x 10~ 3 GeV" 1 to the SM. 

The results in Eq. [11] and Eq. [12] provide a rough estimate for the sensitivity of the 
Tevatron to the CP violating anomalous couplings. The existing Tevatron samples of ti 
events with a double b tag are of the order of 1000 events |l5| and this leads to a 3a 
statistical sensitivity to d t and f t of order 1. To account for background, we notice that: a) 
the experimental cuts to select the ti events are the same that will be used for a CP violation 
study , and b) all the known background processes are CP conserving. The net effect of the 
background (apart from possible systematic errors that must be studied by the experiments) 
is to dilute the asymmetries by a factor (B + S)/S. The numerator in A\ does not get 
additional contributions from the background; but the denominator, which counts the total 
number of events, does. Similarly, the statistical sensitivity decreases by a corresponding 
factor \/ (B + S)/S. For samples with roughly the same number of background (B) and 
signal (S) events this amounts to factors of two. 
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FIG. 3: Differential distributions da/d0 5 for (a) d = 5 x 10 3 GeV 1 and (b) 



/sin, 



5 x 10~ 3 GeV" 1 compared to the SM. 5 is given in units of m\ 



The gluon fusion initiated dilepton channel at y S = 14 TeV leads to A\ = 1.17 dt 
Comparing this number to those in Eq. CCD and in Eq. [12] we see that the dilepton process 
at LHC is an order of magnitude more sensitive than the lepton plus jets or all hadronic 
channels at the Tevatron. This is due to two reasons: first the gluon fusion initiated process 
is more sensitive to the anomalous couplings. For the case of the Tevatron, only about 15% 
of top-quark pairs are produced via this mechanism. Numerically we have seen that if we 
restrict the top-quark pair sample to that originating from gluon fusion, the asymmetries 
increase roughly by factors of three. The second reason is that the di-lepton channel is more 
sensitive to these anomalous couplings. 

In addition, the statistical sensitivity of the Tevatron with 1000 events is about five times 
below that of a 10 fb" 1 run at LHC which would produce about 23000 ti dimuon events 
after the cuts in Eq. |8] and Eq. [9] are applied. Nevertheless, a study with the available 
Tevatron data would be extremely valuable in understanding the role of systematic errors 
in measuring T-odd asymmetries. 

We have performed a series of checks on our numerical analysis as follows. First, we 
evaluate the asymmetries for a few values of the anomalous couplings to check that they 
scale linearly. Second, when the estimated asymmetry is small compared to the 3a statistical 
uncertainty, we repeat the estimate with larger event samples and/or larger values of the 
anomalous coupling to distinguish between zero asymmetries and numerically small ones. 
Third, for ti pair production at the Tevatron the parton process with uu quarks in the initial 
state dominates. We have therefore estimated the asymmetries using this parton process 
only, finding numbers within 10% from the ones obtained when all initial qq and gg states 
are included. 
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IV. SUMMARY AND CONCLUSION 



We have studied the sensitivity of the Tevatron to CP violating anomalous top-quark 
couplings including its chromo-electric dipole moment d. To this effect we have presented 
a numerical implementation of the results in Ref. jsj using Madgraph for event generation 
at the parton level. We have considered processes corresponding to ti events in the lepton 
plus jets and all hadronic channels with two b-tags. In order to generate a statistically 
clean sample we have performed our numerical simulation for a rather large value of the 
anomalous couplings (5 x GeV -1 ). Using the fact that all the asymmetries are linear in 
the anomalous couplings we present our final results as equations in terms of these couplings, 
in Eq. dU H 

Numerically, we find a 3a statistical sensitivity to couplings of order one when normalized 
to the top-quark mass: dt and ft, Eq. [TUJ This sensitivity is about two orders of magnitude 
below what can be accomplished at the LHC with 10 fb" 1 . These results are based on the as- 
sumption that there will be one thousand reconstructed ti events. There could be additional 
inefficiencies in the reconstruction of our specific observables that must be addressed by a 
careful experimental study. A few comments are in order: models available in the literature 
to estimate these anomalous couplings typically yield values too small to be observed at the 
Tevatron; specific models with new sources of CP violation may give contributions to the 
observables we study that cannot be parametrized by the anomalous couplings. 

With a long term goal of searching for CP violation in tt events at the LHC, it is an 
important exercise to analyze the available Tevatron data and we urge our experimental 
colleagues to carry out this study. 
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Appendix A: T-odd correlations 

The spin and color averaged matrix element squared that contains the T-odd correlations 
can be easily obtained from the results in Ref. jsj]. For CP violation in the production process 
qq — > ti they can be written as 

\M\cp = C 1 (s,t,u)0 1 + C 2 (s,t,u)0 2 + C 3 (s,t,u)0 3 , (Al) 
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in terms of the correlations 3 

O x = e(p t ,Pi,PD,PD) 

2 = (t-u)e(p D ,p B ,p,q) 

£>3 = (t-u) (p-p D e(p D ,pt,pt,q)+p-pDe(p D ,pt,Pt,q))- (A2) 

In Eg. lAll and in Eq. IA2l we have used s, t, u, the standard parton level Mandelstam variables 
for tt production. We have also used the sum and difference of parton momenta 

P = P1+P2 

q = Pi ~Pi- (A3) 

Ref. [if explicitly gives the result for the case where both Ws are reconstructed as one 
jet, in which case the form factors are 



Cf(s,t,u) = -^dK bb m t (L-^ +4 <) 
9 V s s J 



16 ~ T ^ m t 



Cl q (s,t,u) = --dK, 



9 """" s 2 

( 

2 



Cf(s,t,u) = S -Cf(s,t,u), (A4) 



and p D = p b , p D = p~ b , 

K bb = (n 2 ay) (2 - <^ 2 (-^ 2 8(pt - ml) 5{p\ - m \). (A5) 

Ref. jsj also indicates how to convert these results into those needed in the case where the 
Ws decay leptonically. For the Tevatron we are interested in two additional cases: 

• Lepton plus jets events. Here one of the W bosons decays leptonically and the other 
one decays hadronically into two jets. The results for a positively charged lepton, £ + , 
follow from the substitutions: 

- In Eq. |A2]pd -> Pe+ and p D -> p d 

- In Eq. |A4] K bb ->■ K u where 



K id = 48 (vr 2 a^ 8 ) (p b ■ p u ) (p b ■ p e ) 



7T \ / 7T 



m t T t J \M W T W 
x 5{p 2 t - tuMpI ' m 2 t )8(p 2 w+ - M 2 w )5{p 2 w . - M 2 W )- (A6) 

For a negatively charged lepton, £~, we need the substitutions: 
— In Eq. IA2l p n —> pg and p^ — > pi- 



Notice that these form factors differ from those defined in Ref. 5] by factors of t — u 
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In Eq. |Al K bb -> K M where 



Km = 48 (vr a s g ) (p 6 • p u ) (p b ■ p v ) 



7T \ / 7T 



m t T t J \M W T W 
x 5(p 2 - m t 2 )<5(p 2 - m 2 )^ - M 2 w )8(p 2 w - - M 2 W )- (A7) 

For the all hadronic decay the following substitutions are required: 

- In Eq. |A2]p D -> pj and ps -> p d 

- In Eq. [53] if 66 ->■ where 



if dd = 144 (ir a s g ) (p 6 ■ p u ) (pg • p € 



7T \ / 7T 



m t T t J \M W T W 
x 5(p 2 - m 2 )5(p 2 - m?)*(& + - M^)5(p^ - M 2 ,); (A8) 

In our numerical implementation we rewrite all delta functions as the respective Breit- 
Wigner distributions behind them, for example: 

KpI - mi) -> . 2 * ■ (A9) 



mJV v ^ l ' (pt~mf)^ + Timi 

When CP violation occurs in the decay vertex, the spin and color averaged matrix element 
squared containing the T-odd correlations was written in Ref. jsj as 4 

\M\t = fsin((f) f + 5f)e(p t ,p b ,p e +,Q t ) + fsin((f) f - 5 f )e(pt,p- b ,p e -,Qt). (A10) 

The Q t and Qi that appear in Eq. IA10I are linear combinations of available momenta 
and act as spin analyzers for the t and i respectively. For tt pairs produced by light qq 
annihilation and when both W's decay leptonically, they are 

Qf = K u {(4sm 2 + (t - u) 2 - s 2 )p E - + 2(sp e - ■ (p t - pi) - (t - u)p e - ■ q)pi 

+ 2((t - u)p e - ■ {p t + Pt) - sp t - ■ q)q] 



Qf = K u {(4sm 2 + {t- u) 2 - s 2 )p £ + - 2(sp t + ■ (p t - pi) - (t - u)p e + ■ q)p t 



16m t 

- 2((t - u)p i+ ■ (p t + pt) + sp e+ ■ q)q} . (All) 
The necessary replacements to obtain the results relevant for us are: 

• Lepton plus jets events with a positively charged lepton, £ + : 

- In Eqs. lATOl \KH\p,- ->■ p d 

- In Eq.EHfe -+K U . 

• For a negatively charged lepton, £~ , 



Note that there is a typo in Ref. |5[ where 4>f and 5f are reversed 
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- In Eq.ElEn]pf+ ->PJ 

- In Eq.\EmK u ^ K M . 

• For the case in which both W's decay into two jets: 

- In Eq. |A10j \All\p e - ->■ p d , p e + ->■ pj 

- In Eq.lMT]^^^. 

Corresponding changes are needed for the gluon fusion processes to the results in Ref. jsl 



[1] See D. Atwood, S. Bar-Shalom, G. Eilam and A. Soni, Phys. Rept. 347, 1 (2001) 
arXiv:fiep-ph/0006032 1; G. Brooijmans et al, larXiv: 0802.3715 [hep-ph]; and references 
therein. 

[2] J. F. Donoghue and G. Valencia, Phys. Rev. Lett. 58, 451 (1987) [Erratum-ibid. 60, 243 
(1988)]. 

[3] J. P. Ma and A. Brandenburg, Z. Phys. C 56, 97 (1992); W. Bernreuther, O. Nachtmann, 
P. Overmann and T. Schroder, Nucl. Phys. B 388, 53 (1992) [Erratum-ibid. B 406, 516 
(1993)]; A. Brandenburg and J. P. Ma, Phys. Lett. B 298, 211 (1993); D. Atwood, A. Aeppli 
and A. Soni, Phys. Rev. Lett. 69, 2754 (1992); W. Bernreuther and A. Brandenburg, Phys. 
Rev. D 49, 4481 (1994) [arXiv:hep-ph/9312210 |; S. Y. Choi, C. S. Kim and J. Lee, Phys. 
Lett. B 415, 67 (1997) [arXiv:hep-ph/9706379] ; H. Y. Zhou, Phys. Rev. D 58, 114002 (1998) 
|arXiv:hep-ph/9805358| ; J. Sjolin, J. Phys. G 29, 543 (2003); J. A. Aguilar-Saavedra, J. Car- 
valho, N. Castro, A. Onofre and F. Veloso, Eur. Phys. J. C 53, 689 (2008) [ arXiv:0705.304T1 
[hep-ph]]; Z. Hioki and K. Ohkuma, larXiv:0910.3049l [hep-ph]; D. Choudhury and P. Saha, 
larXiv:0911.50TBl [hep-ph]. 

[4] D. Atwood and A. Soni, Phys. Rev. D 45, 2405 (1992); D. Chang, W. Y. Keung and I. Phillips, 
Nucl. Phys. B 408, 286 (1993) [Erratum-ibid. B 429, 255 (1994)] [arXiv:hep-ph/9301259] ; 
S. Y. Choi and K. Hagiwara, Phys. Lett. B 359, 369 (1995) |arXiv:hep-ph/9506430| ; P. Poulose 
and S. D. Rindani, Phys. Rev. D 57, 5444 (1998) [Erratum-ibid. D 61, 119902 (2000)] 
[arXiv:hep-ph/9709225] ; S. D. Rindani and M. M. Tung, Eur. Phys. J. C 11, 485 (1999) 
[arXiv:hep-ph/9904319] ; S. D. Rindani, Pramana 54, 791 (2000) [ jarXiv:hep-ph/0002006] ; 
B. Grzadkowski and Z. Hioki, Phys. Lett. B 476, 87 (2000) |[ arXiv:hep-ph/9911505] ; 



14 



Z. H. Lin, T. Han, T. Huang, J. X. Wang and X. Zhang, Phys. Rev. D 65, 014008 (2001) 
[arXiv:hep-ph/0106344] ; S. D. Rindani, Pramana 61, 33 (2003) [arXiv:hep-ph/0304046| . 

[5] O. Antipin and G. Valencia, Phys. Rev. D 79, 013013 (2009) |arXiv:0807.1295 [hep-ph]]. 

[6] S. K. Gupta, A. S. Mete and G. Valencia, Phys. Rev. D 80, 034013 (2009) [arXiv: 0905 .10741 
[hep-ph]]. 

[7] S. Dawson and G. Valencia, Phys. Rev. D 52, 2717 (1995) |arXiv:hep-ph/9504209| ; S. Dawson, 
X. G. He and G. Valencia, Phys. Lett. B 390, 431 (1997) |arXiv:hep-ph/9609523] ; J. Kumar, 
A. Rajaraman and J. D. Wells, Phys. Rev. D 78, 035014 (2008) |arXiv:0801.289T1 [hep-ph]]; 
T. Han and Y. Li. larXiv:0911.2933l [hep-ph], 

[8] P. Langacker, G. Paz, L. T. Wang and I. Yavin, JHEP 0707, 055 (2007) 
[arXiv:hep-ph/0702068] ; A. Szynkman, K. Kiers and D. London, Phys. Rev. D 75, 075009 
(2007) | arXiv:hep-ph/0701165] ; K. Kiers, A. Szynkman and D. London, Phys. Rev. D 74, 
035004 (2006) |arXiv:hep-ph/0605123| . 

[9] M. B. Gavela, F. Iddir, A. Le Yaouanc, L. Oliver, O. Pene and J. C. Raynal, Phys. Rev. D 
39, 1870 (1989); M. P. Kamionkowski, Phys. Rev. D 41, 1672 (1990). 
[10] F. del Aguila and J. A. Aguilar-Saavedra, Phys. Rev. D 67, 014009 (2003) 
|arXiv:hep-ph/020817l1 . 

[11] T. Stelzer and W. F. Long, Comput. Phys. Commun. 81, 357 (1994) [arXiv:hep-ph/9401258] . 
[12] J. Alwall et al, JHEP 0709, 028 (2007) |arXiv:0706.2334l [hep-ph]]. 

[13] J. Alwall, P. Artoisenet, S. de Visscher, C. Duhr, R. Frederix, M. Herquet and O. Mattelaer, 

AIP Conf. Proc. 1078, 84 (2009) |arXiv:0809.24l0l [hep-ph]]. 
[14] G. Mahlon and S. J. Parke, Phys. Rev. D 53, 4886 (1996) [arXiv:hep-ph/9512264| . 
[15] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. D 76, 072009 (2007) |arXiv:0706.3790l [hep- 

ex]]; V. M. Abazov et al. [DO Collaboration], Phys. Lett. B 668, 98 (2008) |arXiv:0804.3664l 

[hep-ex]]. 



15 



